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ABSTRACT 

We present results on the Fundamental Plane of early-type galaxies in the rich cluster 
CL0024+ 16 at z=0.391. The internal velocity dispersions of the galaxies have been 
measured from a 19 hour integration at the MMT. The photometric parameters of 
y—i the galaxies have been measured from deep HST images. The galaxies satisfy a tight 

j Fundamental Plane relation which is similar to that at low redshift. The scatter is 15 % 

C$ in log r e , also very similar to that at low redshift. The data show that massive early-type 

galaxies existed at 2=0.4, and extend earlier studies of the luminosities and colors of 
early-types in rich clusters. 

The evolution of the M / L ratio is derived from the Fundamental Plane, by a com- 
parison with Coma. The M/L ratio increases by 31 ± 12 % between z = 0.391, and 
z = 0.023. The evolution is low when compared to models for stellar populations. The 
(f-) expected evolution depends on the IMF, qo, and the formation redshift of the galaxies. 

\Q The data are in agreement with high formation redshifts. The modeling is still un- 

certain, however, because of various possible biases. The most serious bias may be the 
C*~) progenitor bias: if the progenitors of some current day early-types are spirals at z = 0.4, 

they would not be included in the sample, and the sample would be biased towards the 
oldest galaxies. 

More data are needed to measure the evolution of the Fundamental Plane more 
precisely, and its scatter. There is a hint that the form of the Fundamental Plane 
Cln changes with redshift, and this needs to be determined better. Deeper samples on more 

clusters would be valuable. Studies of the richest nearby clusters may help to test the 
underlying hypothesis that the Fundamental Plane is identical in all clusters. 

^ Key words: galaxies: evolution, galaxies: elliptical, galaxies: mass-to-light ratios, 
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1 INTRODUCTION 

There is substantial evidence that galaxy evolution is a com- 
plex process. In contrast to the earliest ideas, involving an 
early formation and subsequently a smooth and regular dim- 
ming of the stellar light (e.g., Searle, Sargent, & Bagnuolo 
1973), it is now believed that galaxies may evolve in an ir- 
regular way. Interactions, star bursts, and infall of material 
may drive the evolution of a typical galaxy. Observational 

* Observations reported here are partly obtained at the Multiple 
Mirror Telescope, a joint facility of the Smithsonian Institution 
and the University of Arizona. 

t Visiting Astronomer, Kitt Peak Observatory, a division of 
National Optical Astronomy Observatories, which is operat- 
ed by AURA, Inc., under contract with the National Science 
Foundation. 



evidence comes from the evolution of cluster galaxies (e.g., 
Butcher & Oemler 1978, 1984, Dressier & Gunn 1983, Couch 
& Sharpies 1987), and studies of field galaxies (e.g., Broad- 
hurst, Ellis, & Shanks 1988, Lilly et al. 1995). Similarly, hi- 
erarchical formation scenarios predict that galaxy evolution 
may be an irregular process (e.g., White & Frenk 1991). 

Such complex evolution makes the interpretation of the 
observations more ambiguous, as the luminosities of galaxies 
can vary rapidly due to star bursts, and the morphologies of 
galaxies can change with time. As a result, the progenitors 
of current galaxies may have had different morphologies and 
masses. 

Detailed studies of the evolution of galaxy morphology, 
luminosity, and mass are necessary to determine the rele- 
vance of these processes. Such work has now become possi- 
ble with the high resolution imaging capability of the Hubble 
Space Telescope (e.g., Dressier et al. 1994), and deep spec- 
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troscopy from the ground. An essential aspect of such stud- 
ies is the measurement of the evolution of the masses and 
mass-to-light ratios (M/L) of galaxies. These measurements 
constrain the mass evolution of galaxies, and thereby infall 
and merging. They are necessary for a correct interpretation 
of the evolution of the luminosity function. 

As the total masses of galaxies are very difficult to mea- 
sure, relations such as the Tully-Fisher relation (Tully & 
Fisher 1977), the Faber- Jackson relation (Faber & Jackson 
1976), or the Fundamental Plane (Dressier et al. 1987; D- 
jorgovski & Davis 1987) can be used to constrain the mass 
evolution of galaxies. The small intrinsic scatter of these re- 
lations makes them very useful for studies of evolutionary 
effects. 

Here we present results on the Fundamental Plane re- 
lation at high redshift. The Fundamental Plane is a re- 
lation between effective radius r e , effective surface bright- 
ness I e , and central velocity dispersion a of the form r e cx 
0.1.24 j-o.82 (Bender, Burstein, & Faber 1992, J0rgensen, 
Franx & Kjaergaard 1995c [JFK95c]). As shown originally 
by Faber et al. (1987), the implication of the relation is that 
the M/L ratio of galaxies is well-behaved, and of the form 
M/L oc r° 22 (7° 49 oc M° 24 . Recent studies indicate that the 
scatter in the M/L for cluster galaxies is small (23%, e.g., 
Lucey et al. 1991, JFK95c). The small scatter makes the 
Fundamental Plane especially useful to measure the evolu- 
tion of galaxies with redshift, as small samples can provide 
useful answers. Furthermore, biases due to sample selection 
effects are small, if the scatter is small. The study of the 
Fundamental Plane at higher redshift has recently become 
possible with high throughput, high spectral resolution spec- 
trographs (e.g., Franx 1993a, b). 

In this paper, we report on measurements of velocity 
dispersions and structural parameters of galaxies in the rich 
and concentrated cluster CL 0024 + 1654 at z = 0.39. The 
Fundamental Plane is derived, and the implications for the 
evolution of the mass-to-light ratio are discussed. Sections 2 
and 3 present the spectroscopy and the photometry. Those 
readers who are not interested in the details of the data 
analysis can proceed directly to sections 4 and 5, where the 
Fundamental Plane is presented and discussed. 



2 SPECTROSCOPY 

Objects in the CL 0024 + 1654 field were observed with the 
Multiple Mirror Telescope (MMT) in September 1992. A 
mask with 14 slitlets was used. Each slitlet had a width of 
1.6 arcsec, and variable length. Objects were selected on the 
basis of their R band flux (restframe ~ B). The slit posi- 
tions were optimized to include as many of the brightest 20 
galaxies as possible, and fainter galaxies were used to fill in 
remaining positions. The total integration time was 70 200 s, 
consisting of thirteen 5 400 s exposures in three nights. One 
object turned out to be a star, the thirteen galaxies are listed 
in Table 1, along with descriptions of low resolution spectra 
taken by Dressier, Gunn, & Schneider (1985) [DGS]. In this 
paper, we will use the numbering system of DGS through- 
out. 



2.1 Reduction 

Each slitlet was treated as a separate, long slit spectrum. For 
the reduction, we primarily used the standard IRAF tasks. 
Some tasks were written for the removal of the cosmic ray 
hits. 

After bias subtraction a smooth gradient was present 
in the residual bias frames. The peak-to-peak variation in 
this gradient proved to be small (0.4%), and stable in time. 
We measured the gradient by fitting a low-order function to 
the residual of an average bias frame of each night, and we 
subtracted this fit from the data. The remaining variations 
in the individual bias frames are due to noise. Internal flat- 
fields were used to correct for the pixel- to- pixel variations in 
the CCD response. The position of the slits on the chip was 
relatively stable during the observations: total shifts were of 
order ~ 2 pixels during each night. Several averaged internal 
flatfields were created for each night to minimize the effect- 
s of these shifts. Since internal flatfields were taken before 
and after each 5 400 s exposure, small changes in the flat- 
fields over time (e.g., the accumulation of dust grains) could 
be monitored well. 

Cosmic ray hits were removed in several steps. First, 
exposures in each night were median filtered. The median 
frames were subtracted from the exposures, and the resid- 
ual galaxy light and the residuals of skylines were modeled 
with appropriate functions. Cosmic ray hits were identified 
automatically as significant deviations from the function fit- 
s. The results from the program were inspected by eye, and 
corrected if necessary. Left over cosmic rays (~ 10 per expo- 
sure) and bad pixels were removed interactively. Finally, the 
separate exposures of each night were aligned and combined, 
with optimal weighting. 

For each night, a wavelength solution was obtained from 
exposures of He-Ne-Ar lamps. The frames were transformed 
to a common log A scale. The He-Ne-Ar exposures were re- 
duced and combined in the same way as the galaxy expo- 
sures. The response function along each slitlet was deter- 
mined from twilight sky exposures. After removal of the re- 
sponse function systematic variations in the spatial direction 
were < 1 %. The sky spectrum was derived from the spec- 
tra at the edges of each slitlet. After sky-subtraction, the 
average frames of the three nights were combined to form 
averaged, wavelength calibrated spectra of the 14 objects. 

2.2 Instrumental Resolution 

The usual procedure to determine velocity dispersions is 
to compare a galaxy spectrum with a star spectrum tak- 
en through the same spectrograph with the same setup. 
This ensures that the instrumental resolution of the star 
spectrum is the same as the instrumental resolution of the 
galaxy spectrum. Spectroscopy of high redshift galaxies does 
not allow this method, as the star and galaxy spectra can- 
not be obtained with the same setup. For the observation 
of stellar spectra, the grating angle needs to be changed to 
cover the same rest wavelengths as the spectra of the high 
redshift galaxies. Since the resolution is roughly constant in 
A, it will change significantly in kms -1 (Franx 1993a, b). 
Furthermore, the resolution is usually dependent on the po- 
sition on the CCD. Hence an absolute determination of the 
resolution, over the whole CCD, is necessary. 
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Table 1. Spectroscopic Parameters 



DGS a 


BO b 


PK C 


Spectral Description 1 ^ 




5(7 


R>emarks 


84 


15 


122 


Weak H,K,G,Bal abs 


- 


- 


Insufficient S/N 


111 


26 


105 


E-type 


156 


28 




130 




68 




243 


22 




137 


42 


80 








Insufficient S/N 


158 


62 


87 


E-type 


317 


23 




161 


65 


86 


E-type 


275 


20 




162 


72 


108 


E-type 


167 


19 




169 


71 


88 


E-type 


342 


25 




186 


87 


61 


E-type: 


382 


34 




194 


96 


85 


O II, Bal abs 






Insufficient S/N 


202 


100 


70 


Weak E-type 


248 


20 




218 








145 


33 




229 


125 


64 








Insufficient S/N 



Notes: 

a Dressier, Gum & Schneider (1985) [DGS] number. 

6 Butcher & Oemler (1978) number. 

c Pickles & van der Kruit (1991) number. 

d Description of low resolution spectra from DGS. 



We determined the instrumental resolution from the He- 
Ne-Ar exposures, by fitting Gaussians to the emission lines. 
The fits were made for each slitlet, and for each observ- 
ing night separately. Then, smooth functions were fit to the 
measured dispersions. The procedure is illustrated in Fig. 1. 

The instrumental resolution proved to be approximate- 
ly constant at <Ti nstr = 1.4 A (70kms _1 at 6000 A) for 
A ;> 6100 A. At shorter wavelengths, the instrumental res- 
olution is larger: it rises to ~ 150 km s -1 at 5500 A. The 
resolution also depends on the spatial coordinate: for the 
top and bottom slitlets cr; ns t r is ~ 10 % larger than for the 
middle slitlet. The resolution was stable from night to night. 

Since each cluster exposure was bracketed by He-Ne- 
Ar exposures, the resolution could be determined through 
the nights. The resolution proved to be stable during each 
night, and consistent with the results derived from the mean 
He-Ne-Ar exposure. 

The profile of a single line is determined by the instru- 
mental characteristics, and can deviate from a Gaussian. 
This is especially likely to occur when the camera is out of 
focus on parts of the CCD. We verified the camera focus, and 
the shapes of the He-Ne-Ar lines by means of two exposures 
taken through two Hartman masks. These masks mask out 
either half of the parallel beam in the spectrograph. If the 
camera is out of focus, then there will be an offset between 
the two exposures in the wavelength direction. The offsets 
were small (~ 0.6 A) in the spectral region that we used. 
Furthermore, we modeled the total instrumental resolution 
by taking the sum of the two separate Hartman exposures. 
The He-Ne-Ar lines on each Hartman exposure were fitted 
separately, and the instrumental resolution was now defined 
as the sum of two Gaussians, with different widths, and an 
offset. The resulting line profiles were very similar to those 
derived from a single Gauss fit, and the profiles did not af- 
fect the final velocity dispersions of the galaxies. This gives 
us confidence that the modeling is adequate. 

The instrumental resolution can also be determined 
from sky emission lines, present in the galaxy spectra. This 



method has the advantage that the actual galaxy exposures 
are used, but suffers from the fact that many sky lines are 
blends. Gaussians were fit to bright lines which appeared to 
be no blends, using the final, summed spectra. No system- 
atic differences between the thus determined resolution and 
the results from the He-Ne-Ar exposures were found. Since 
all determinations of the instrumental resolution give con- 
sistent values, we are confident the instrumental resolution 
is correctly modeled. 

High spectral resolution template stars were kindly 
made available by R. van der Marel. The observations were 
done May 22, 1992, with the William Herschell Telescope 
(WHT) on La Palma, with the blue arm of ISIS (van der 
Marel et al. 1994). The spectral types of the 16 stars range 
from GO to M0. The spectral resolution was measured again 
from comparison lamp spectra, similar to our measurements 
of the resolution of the MMT data. The spectral resolution 
is constant as a function of wavelength, but depends on the 
width of the slit used in the observations: a* = 0.35 A (l" 
slit) and a* = 0.62 A (l'.'8 slit). The spectra cover the spec- 
tral range 4215 - 4615 A. 

The spectra were smoothed, rebinned, and smoothed 
again to have identical resolution as the galaxy exposures. 
Since the resolution of the galaxy spectra depends on the 
position on the chip, a separate template star was produced 
for each galaxy spectrum. 

We checked our method for consistency by comparing 
velocity dispersions obtained from template star spectra ob- 
served through the wide and the narrow slit. The differences 
in the output dispersion are of order 1 %, and are not sys- 
tematic. Hence our procedure correctly compensates for the 
differences in spectral resolution. 

Eight additional template stars were made available by 
D. Fisher. They were observed May 16, 1993, with the Kast 
Spectrograph on the Lick Observatory 3 m telescope. The 
Lick spectra have a lower spectral resolution (a* = 1.36 A), 
but cover a larger spectral range (4215 - 5617 A) than the 
WHT stellar spectra. The stars were redshifted and rebinned 
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Figure 1. The instrumental resolution as a function of wave- 
length and position on the CCD. The instrumental resolution has 
been measured from He-Ne-Ar exposures. Solid symbols indicate 
the resolution for the central slitlet, open symbols the resolution 
for the slitlet closest to the edge of the chip. The deviant point 
at 6050 A is a doublet. The resolution is highest in the centre of 
the chip. The lines show the smooth functions that were fit to 
the measurements. We used the wavelength interval from 5895 to 
6400 A to measure the velocity dispersions. 

to mimic the MMT galaxy data. Since the instrumental res- 
olution is comparable to the instrumental resolution of our 
data, no further transformations were performed. 

2.3 Fourier Fitting 

Velocities and velocity dispersions were determined with 
the Fourier Fitting method (Franx, Illingworth, & Heckman 
1989 [FIH]). This method minimizes the residuals of the fit- 
ted stellar spectrum to the galaxy spectrum 

N 

x i = J2(G,-[BoS],f/s 2 (1) 

8 = 1 

where G is the galaxy spectrum, B is the broadening func- 
tion, S is the stellar spectrum, o denotes a convolution, and s 
is the noise in the galaxy spectrum. The fit is done in Fouri- 
er space, which speeds up the computation. The technique 
is described in detail in FIH. 

Several other methods exist for the determination of 
velocities and velocity dispersions from galaxy spectra (see, 
e.g., Rix & White 1992). One of our objectives is to com- 
pare our results with a low redshift comparison sample. 
J0rgensen, Franx, & Kjaergaard (1995b) [JFK95b] have used 
the Fourier Fitting method for the determination of velocity 
dispersions of a large sample of galaxies, and they found re- 
sults consistent with earlier literature data. We use the same 



program in a very similar setup, to limit any systematic er- 
rors. 

The spectra were summed in the spatial direction to in- 
crease the S/N ratio. The three central spectra were added to 
give an effective aperture size of l'.'8 x l'.'6. The S/N decreas- 
es very rapidly at larger distances from the galaxy centers, 
and the effective S/N ratio cannot be improved by adding 
more spectra. No weighting was applied in the summation. 

The Fourier program subtracted the continuum of the 
spectra, tapered them at the ends, and calculated the Fouri- 
er transform. The spectra were filtered with a low bandpass 
filter, and the fit was performed. The filter was set to give 
zero power at Fourier wavenumbers zero to l , and to rise s- 
lowly to unity at 2l . We chose l a equivalent to 110 A -1 . This 
is identical to the value used by JFK95b when expressed in 
units of inverse km s -1 . 

Since the spectra are not flux calibrated, the shape of 
the continuum is partly determined by the response function 
of the detector. The continuum shapes of the galaxy spectra 
and the star spectra are somewhat different, because they 
were taken with different instruments. The error introduced 
by this difference was estimated by dividing out the slope 
of the stellar spectrum, and comparing the results from the 
fitting program. The resulting rms differences in a are ~ 1 %. 
Hence the fitting technique is insensitive to the shape of the 
continuum. 

Various wavelength intervals for the fits were tried. The 
blue limit of the (redshifted) WHT and Lick star spectra is 
5863 A. Inclusion of the interval 5863 - 5895 A resulted in 
a large increase in \ 2 f° r several spectra, due to the pres- 
ence of a strong sky emission line at 5890 A. Therefore, only 
the part of the spectrum redward of 5895 A (rest wavelength 
4238 A) was used in the fits. For the fits with the WHT tem- 
plate stars, the upper limit of the wavelength interval was 
determined by the red limit of the redshifted star spectrum 
at 6420 A. When fitting the WHT stars, we thus used the 
interval 5895 - 6400 A. The Lick stellar spectra cover a larg- 
er wavelength range. The red limit is not determined by the 
stellar spectrum, but by the galaxy spectrum. This limit is 
different for each galaxy; it is determined by the layout of 
the multi-aperture mask. Typically it is ~ 6530 A. The re- 
sults of the Fourier Fitting did not change when a red limit 
of 6400 A was used. 

The match between the template star spectrum and the 
galaxy spectrum will never be perfect, since the galaxy spec- 
trum consists of a mix of stellar populations. In practice, 
template mismatch can be minimized by fitting template s- 
tars of various spectral types and using a star of the spectral 
type that gives the lowest \ 2 residual. 

Since the total number of template stars is 24, and the 
number of galaxies is 13, the best fitting templates could 
be determined reliably. The \ 2 residual has a minimum at 
spectral types G9 - Kl, and it rises gradually for cooler 
or hotter stars. The dependence of \ 2 on spectral type is 
strongest for the galaxies with the highest S/N. 

Three of the WHT stars and five of the Lick stars have 
spectral types in the range G9 - Kl. The consistency of the 
results of the fitting program was estimated by comparing a 
and 8a obtained with these template stars. The agreemen- 
t between the Lick stars and the WHT stars is good. For 
the highest S/N spectra the distributions of the individual 
measurements of a and 8a are indistinguishable; for a few 
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galaxies the differences are up to 10 %. There is no system- 
atic effect, i.e., the differences average out for the galaxies. 
The agreement gives confidence in the stability of the re- 
sults. Given the fact that the instrumental resolution of the 
Lick stars is not precisely that of the galaxy spectra (espe- 
cially blueward of 6100 A, where the resolution of the galaxy 
spectra is no longer constant), the good agreement was not 
expected a priori. 

Since the spectral resolution of the WHT templates is 
the same as the instrumental resolution of the galaxy spectra 
over the whole wavelength region used in the fits, we chose 
a good fitting WHT template (the K0 III giant HD 132737) 
for the final determinations of the velocities and velocity dis- 
persions. Several galaxy spectra did not have sufficient S/N 
to derive a proper velocity dispersion. These galaxies were 
omitted from the list, as low S/N spectra produce biased 
estimates of the velocity dispersion (see, e.g., FIH). Fig. 2 
shows the spectra of the galaxies for which reliable velocity 
dispersions could be determined. The position of the bright- 
est sky emission line is indicated by bars. 

Early-type galaxies have radial gradients in the velocity 
dispersion (e.g., Illingworth 1981, Davies et al. 1983, FIH, 
Jedrzejewski & Schechter 1989, van der Marel, Binney, & 
Davies 1990, Fisher, Illingworth, & Franx 1995). Therefore, 
the measured velocity dispersion through an aperture de- 
pends on the distance of the object and the aperture size. 
For a meaningful comparison between galaxies, the mea- 
sured velocity dispersions must be corrected to a common 
physical aperture size. JFK95b have established aperture 
corrections from empirical models, based on published pho- 
tometry and spectroscopy of 51 galaxies. They found that 
the correction to a physical aperture size of d covv kpc can be 
approximated well by a power law: 

log <T cor r = log <T ob s + 0.04 (log d obs - log (2 C orr), (2) 

where d b s is the physical size in kpc of the aperture used in 
the observations. The correction is a simple power-law, and 
therefore the effective radius does not enter the expression 
if the aperture is corrected to a fixed physical size in kpc. 
We use the same normalisation as JFK95b: a 3'.'4 diameter 
circular aperture projected on a galaxy in Coma. In our 
observations, the effective aperture was l'.'6 x l'.'8. This is 
equivalent to a circular aperture of l'.'96 diameter (JFK95b). 
Using Eq. 2, we find that the correction to the distance of 
Coma and an aperture size of3'.'4is <r corr = 1.072 <r bs, where 
we used a value of go = 0.5. 

Corrected velocity dispersions are listed in Table 1. The 
errors are the formal errors from the Fourier fitting program. 
Systematic errors in the velocity dispersion may be produced 
by two effects: spectral type mismatch, and incorrect filter- 
ing. We estimated these from the results obtained with tem- 
plate stars with different spectral type, and by changing the 
filter l by a factor of 1.5. Each effect is on the order of 3%, 
and we estimated the combined systematic error at 5%. 



3 PHOTOMETRY 

Hubble Space Telescope (HST) Wide Field Camera (WFC) 
images of CL 0024 + 1654 were obtained from the Space 
Telescope Science Institute (STScI) archive. The cluster was 
observed through the F702W filter on two occasions: Octo- 



ber 6, 1991 and October 1, 1992. The pointings were offset by 
62". Therefore, approximately one half of the objects were 
imaged in both observations. Each observation consists of six 
separate exposures. The images of the 1992 observation are 
not aligned: between the third and fourth exposure the tele- 
scope moved 3'.'7. Total exposure times are 11 200 seconds 
(1991 observation) and 13 100 seconds (1992 observation). 

3.1 Reduction 

The pipeline processing was performed at STScI. For each 
observation we checked whether the recommended calibra- 
tion files were used. The pipeline bias frames were different 
from the recommended frames. We retrieved the used and 
recommended bias frames from the archive and subtracted 
them; no differences were found. Bad columns were interpo- 
lated. Cosmic rays were removed in three steps. In the first 
step, individual exposures of the same field were compared, 
and pixel values which were significantly different from the 
mean were rejected. Then, the exposures were averaged, and 
left-over hits were removed using the COSMICRAYS task in I- 
RAF. Finally, the images were carefully checked by eye for 
both left over cosmic rays and pixels which were wrongly 
recognized as cosmic ray hits (e.g., the central pixels of s- 
tars are sometimes marked as cosmic rays by the software). 

The resulting images show a large (~ 20 %) gradient 
in the background across the four chips. This is probably 
caused by non-uniformity in the throughput of the F122M 
filter, which was used as a neutral density filter when Earth 
flatfields were taken (Faber 1992). Furthermore, "doughnut" 
like features are present (~ 5 % variations), probably images 
of the secondary mirror assembly (see, e.g., Ratnatunga et 
al. 1994). 

First, we tried to remove the gradient by applying a cor- 
rection flatfield, constructed from the "superskyflats" in the 
F555W and F785LP filters created by the Medium Deep 
Survey (MDS) team (Phillips et al. 1994; Ratnatunga et al. 
1994): the MDS flatfields were divided by the stacked earth 
flatfield in the same filter, yielding a F122M filter flatfield. 
These flatfields removed the "doughnuts" , but did not re- 
move the gradient satisfactorily. 

Therefore, we created correction flats from the expo- 
sures themselves, by fitting a low order function to the back- 
ground, carefully avoiding all objects. The correction frames 
were normalized to the center of WF2. These flats do not re- 
move the "doughnuts" , but do reduce the large scale fluctu- 
ations to ~ 1.5 %. The doughnuts were modeled and divided 
out separately, using the constructed FY22M flats. 

As an independent test, we obtained images of the clus- 
ter Abell 370 from the STScI archive, taken November 11, 
1992, also through the F702W filter. In these observations, 
the same background structures are present, confirming that 
the source of the gradient is a natfielding problem. After di- 
vision of the Abell 370 images by our correction flats, the 
background was flat within 2%. 

Finally, we combined the two averaged images from the 
1992 observations, which were offset by 3'.'7. Since the shape 
of the HST PSF is dependent on the position on the chip, 
this step requires justification. We have tested the effect on 
galaxy parameters by refitting all (43) objects on chip 1 
of the 1991 observation with PSFs offset by 18 pixels, and 
comparing the resulting parameters with the original ones 
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Figure 2. Spectra of galaxies with sufficient signal to determine internal velocity dispersions. The numbering system is from Dressier, 
Gunn &l Schneider (1985) [DGS]. The lower left panel shows the template spectrum which was used to determine velocity dispersions. 
The units on the vertical axes are instrumental counts. Redshifts and velocity dispersions are listed in Table 1. Prominent spectral 
features are, e.g., the G-band at 4300 A and the lines at Fe 4383, 4531, and Ca 4455 A. 

shifted with the redshift of the galaxies. The F702W band 
for the HST observations is fairly close to the redshifted V 
band, which we denote by "T4". We use a small color correc- 
tion to bring the HST observations on the V z system. The 
steps are outlined in detail in Appendix A. 

In summary, we have used groundbased data to deter- 
mine a proper zeropoint for the HST instrumental magni- 
tudes, and to determine the color correction. We find 

V z = m_F702W-(0.47±0.05)(i?-/-0.82)-(2.04±0.01),(3) 

The uncertainties in the steps are small. The systematic ze- 
ropoint uncertainty in the V z magnitude is 0.03. The noise 
added by the color term is small, on the order of 0.02 mag. 
As a result, the method produces the true V z mag- 



(cf. Sect. 3.3). The average difference in r e is 0.2 ± 0.9% 
and in I e 0.3 ± 1.7%. The difference in the product r e I°' 8 
is < 0.1 %. The differences in galaxy parameters do not cor- 
relate with distance from the center of the chip. Since the 
differences are small and are not systematic, no significant 
error was introduced by combining the six exposures of the 
1992 observation. 



3.2 Zeropoints 

A meaningful comparison between magnitudes of galaxies at 
large z and magnitudes of nearby galaxies requires a trans- 
formation to a common photometric band in the restframe 
of the galaxies. For this purpose, the filterband needs to be 
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nitudes for galaxies. It is very different from applying 
a /('-correction to magnitudes taken in a fixed filter in 
the observers frame. As noted by many authors, such K- 
corrections produce errors due to color variations between 
ellipticals, and because of color evolution of galaxies with 
redshift. As an example, the scatter in V — R for Coma el- 
lipticals is ~ 0.06 magnitudes. Furthermore, the models of 
Worthey (1994) predict an evolution in V - R of ~ 0.10 
magnitudes between z = 0.4 and the present. Therefore, the 
various /('-corrections in the literature can easily produce er- 
rors of order 0.1 mag. Our resulting magnitudes behave like 
the standard /('-corrected magnitudes, i.e., as fluxes, and not 
as flux densities. This is relevant for the surface brightness 
dimming correction. 

We determined the galactic extinction in the R and / 
band from Burstein & Heiles (1984), and Seaton (1979). 
We obtain an extinction of 0.09 in R, and 0.07 in /. All 
photometry is corrected for this extinction. 

3.3 Analysis 

3.3.1 Fitting Program 

For an accurate measurement of length scales and central 
surface brightnesses it is necessary to account for the Point 
Spread Function (PSF). Especially for the HST, with its 
extended and asymmetrical PSF, an understanding of the 
effects of the PSF is crucial for a reliable determination of 
galaxy parameters. 

We fitted galaxies with a code which uses the full 2D 
information, taking the PSF into account. The program con- 
volves a 2D model of a galaxy with the image of the PSF, 
and iteratively minimizes the residual 

N 

X 2 = ^»,(I,-[Moi'],) 2 , (4) 

8 = 1 

where / is the observed intensity, M is the model, P is the 
PSF, and o denotes convolution. The summation is over all 
pixels within some fixed radius. A weighting function w is 
applied: 



This weighting ensures that the weight of pixels at the effec- 
tive radius is roughly half of that of the central pixels. As a 
stability check, we also fitted objects with a uniform weight- 
ing (i.e., w t = I) and investigated the change in derived pa- 
rameters. The effect of the applied weighting function on r e 
and / e is < 2 %, and the effect on the FP parameter r e I° 8 
< 0.5%. This shows that the derived parameters do not 
depend critically on the applied weighting function. 

The models are parameterized by 7 parameters: x and y 
centre position, effective radius r e , surface brightness at an 
effective radius I e , ellipticity e, position angle, and sky val- 
ue. An additional parameter is object type: exponential, de 
Vaucouleurs r", exponential + r?, or a ^-function, for the 
modeling of stars. All objects are fit simultaneously. For each 
object, a different PSF can be used. The advantage of fit- 
ting to the original data rather than to restored data (e.g., 
CLEANed data) is that the method requires a convolution 
instead of a deconvolution. This ensures that our method 



is repeatable and stable, and gives direct error estimates. 
Furthermore, the correlations between the parameters fol- 
low directly from the least squares fit. This is relevant here, 
because r e and I e are generally strongly coupled in these 
determinations. Fortunately, the correlation is usually par- 
allel to the line of constant r e I°' 8 , and produces very little 
scatter in the Fundamental Plane (e.g., J0rgensen, Franx, & 
Kjaergaard 1993). The other advantage of this procedure is 
that overlapping galaxies are treated properly. 



3.3.2 Results from the Fits 

All objects in the observed CL 0024 + 1654 field with instru- 
mental magnitudes rriF702W < 24.6 were fit. We used Tiny 
Tim 4.0 (Krist 1994) to calculate PSFs with a total size of 
6" x 6". These were calculated for all object positions, to 
account for the dependence of the shape of the PSF on the 
position on the chip. Each galaxy was fit three times: in the 
first two fits, an exponential profile and anr< profile were fit 
respectively, to determine which profile type has the lowest 
X 2 residual. Then, a third fit was performed with the best 
profile type for each galaxy. 

For eight galaxies in our sample an r" profile gives the 
lowest x 2 residual. Galaxy DGS 162 is best fit by an ex- 
ponential profile. We tested how well we can discriminate 
between disk-dominated and bulge-dominated galaxies by 
simulations. We constructed model galaxies, with S/N and 
structural parameters typical for the observations. The mod- 
els cover a range in inclination and bulge-to-disk ratio. We 
found that galaxies with bulge to disk ratios higher than 2/3 
are best fit by r 1//4 law profiles. Hence, galaxies that are best 
fit by an r * profile may have significant disks, contributing 
up to ~ 60% of the light. The S/N is very high for the 
observations, and does not influence this result. 

Fig. 3 (Plate @@) shows the images of the galaxies with 
measured velocity dispersions (cf. Sect. 2.3). For each galaxy, 
the panels show the original image, the residual from the fits 
to the data, and a CLEAN restoration (Hogbom 1974) of the 
image respectively. The CLEANing was performed with a 1.5a 
stopping criterion and a loop gain of 0.02. The results are 
stable with respect to changes in these parameters. Table 2 
lists the values of % 2 , and the effective parameters r e and /t e 
for the galaxies. The value of % 2 is significantly higher than 
I for all galaxies. This is no surprise, as ellipticals are well 
known to have radial variations in ellipticity and position 
angle, and intensity profiles which deviate from an r 1//4 law. 
The galaxies with the largest % 2 are interacting galaxies (see 
3.3.3). 

As an independent check on our procedure, surface 
brightness (SB) profiles were extracted from the original and 
CLEANed images by integrating along ellipses with elliptici- 
ties and position angles as found by our fitting program. The 
profiles are shown in Fig. 4. The dotted lines indicate the 
result from the fit to an r< model; for DGS 162 the dashed 
line indicates the result from the fit to an exponential model. 

We also verified that the direct fitting method did not 
miss flux within an r e . To that end, we derived the fluxes of 
the convolved model and that of the data within an effective 
radius of each galaxy. The differences were small, smaller 
than 0.03 magnitudes. 
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Figure 3. Hubble Space Telescope Wide Field Camera images of galaxies with measured velocity dispersions. Each panel is 13" X 13". 
North is up, East is to the left. For each galaxy, the original data, the residual of the model fits to the data, and a CLEAN restoration is 
shown. Apart from DGS 186, which has a triple nucleus, the galaxies are fit well. The CLEANing tends to sharpen and emphasize features 
already marginally visible in the original data. 
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Figure 4. Surface brightness profiles of galaxies with measured velocity dispersions. Open circles represent the original data, filled circles 
a CLEAN restoration. Dotted lines indicate r 4 profiles as found by our fit to the original data. Some galaxies exhibit deviations from 
an r 4 profile, due to the presence of disks, or other galaxies. Galaxy DGS f 62 is best fit by an exponential profile, indicated with the 
dashed line. Galaxy DGS f 86 has a triple nucleus. 



3.3.3 Morphologies of Galaxies with Measured Velocity 
Dispersions 

Interestingly, several of the SB profiles depicted in Fig. 4 
show significant deviations from a perfect r" law. The pro- 
files are measured along ellipses with constant ellipticity, po- 



sition angle and centre position. The profiles do not change 
when the shapes and the centers of the ellipses are allowed 
to vary. 

We briefly discuss the morphological characteristics of 
the individual galaxies, on the basis of the ellipse fits, the 
CLEAN restorations of the images, and the residuals from 
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the fits to the original images. We note that the distinction 
between ellipticals and SO's is difficult to make for lower 
luminosity galaxies, as face-on disks can be difficult to detect 
(e.g., Rix & White 1990, J0rgensen & Franx 1994). 

DGS 111 - This galaxy has the appearance of an SO close 
to face-on. The inner part of the galaxy is well fit by an r 7 
profile. At r = l'.'S a small companion galaxy causes a bump 
in the profile. The profile in the outer parts (r > 2") is best 
fit by an exponential profile, indicating this galaxy probably 
has a disk. 

DGS 130- The CLEANed image, and the residual map shows 
a tidal tail to the North. The profile changes abruptly at 
r PS l'.'5; the fit represents the inner part well. The outer 
part, containing most of the light, is fit badly. This is where 
the galaxy has an asymmetric extension to the North. We 
conclude that this galaxy is most likely interacting. 

DGS 158- Elliptical. This is one of the four bright galaxies 
in the centre of the cluster. The light profile is consistent 
with an r* law. 

DGS 161 - Elliptical. This is the central galaxy of the central 
group in the cluster. This group consists of four galaxies, 
including DGS 169 and DGS 158. The profile outside of 5" 
is affected by other galaxies. The inner 5" are well fit by an 
r* law. 

DGS 162 - This is an edge-on disk galaxy. The residual 
image and the CLEAN image indicate dust absorption near 
the center and in the disk. The morphological type is SO/S. 
The profile is best fit by an exponential light distribution. 
The galaxy is red (R — I = 0.78); it resembles the bright, red 
spirals observed in nearby rich clusters. Since we intend to 
limit the biases as much as possible, this galaxy was treated 
as an elliptical (i.e., an r< law was fit to the profile). In 
Fig. 4, this fit is indicated with the dotted line. 

DGS 169- Elliptical. The galaxy image is close to the edge of 
the chip, therefore, points at r > 2" have large uncertainties. 
The inner part of the profile is well fit with an r" profile. 

DGS 186- The CLEAN image of DGS 186 clearly shows the 
presence of three nuclei (cf. Fig 3), separated ~ 0'.'5 (3kpc), 
a characteristic of nearby Brightest Cluster Galaxies (e.g., 
Hoessel 1980). The galaxy is not located at the center of the 
cluster, but at a projected distance of 35" (0.22 Mpc). The 
bright nuclei severely affect the fit to the galaxy light profile. 
Therefore, we do not include this galaxy in our analysis. 

DGS 202- This galaxy appears to be an elliptical, or a close 
to face-on SO. It is difficult to distinguish faint ellipticals 
from face-on SO galaxies. The light profile is consistent with 
an r" law. 

DGS 218 - Probably a barred SO galaxy. The galaxy has 
a high ellipticity at r PS 0'.'7. There is a break in the light 
profile at r = 0'.'7, and a position angle twist. The profile is 
exponential in the outer parts, consistent with the presence 
of a disk. 

From the analysis of the HST images, it is clear that var- 
ious types of galaxies are represented in our sample, ranging 
from a disk dominated spiral galaxy (DGS 162) to apparent- 
ly 'bona fide' giant ellipticals like DGS 158. Three galaxies 
show evidence for disks, one galaxy is probably interacting. 
DGS 186 has three nuclei, and is not considered in the anal- 
ysis of the Fundamental Plane. 



Table 2. Photometric Parameters 



Iden a m r Type \ 2 logr e /i^ R c ~ ^ 

n 



111 


19.51 


SO 


4.1 


0.295 


23.28 


24.11 


0.62 


130 


19.47 


Interact 


44 


-0.048 


21.49 


22.37 


0.84 


158 


18.48 


E 


12 


0.773 


23.81 


24.69 


0.83 


161 


18.37 


E 


13 


0.535 


23.21 


24.11 


0.85 


162 


19.48 


SO/S 


12 


0.594 


23.94 


24.81 


0.78 


169 


18.57 


E 


15 


0.348 


22.56 


23.45 


0.86 


186 


19.20 


Triple 


205 


0.601 


23.52 


24.40 


0.81 


202 


19.08 


E/SO 


14 


0.759 


24.08 


24.96 


0.82 


218 


20.55 


SBO 


5.6 


-0.089 


22.06 


22.95 


0.91 



Notes: 

a Gunn r magnitude from Schneider, Dressier & Gunn (1986). 



3.3.4 Internal Accuracy of the Fit Results 

The fact that no systematic signal is present in the residual 
maps (see Fig. 3, Plate @@ ) is reassuring. Also, from Fig. 5 
it is clear that the fit results generally are consistent with 
a CLEAN restoration. Differences between fits and data can 
be attributed to characteristics of the galaxies (e.g., disks 
or tidal features), and do not stem from systematic errors 
in the fitting procedure. Although the formal error in the 
derived parameters is small (typically ~ 1 — 2 %), we suspect 
the true uncertainty is larger, and primarily due to lack of 
knowledge of the HST PSF. 

The shape of the HST PSF is position dependent due 
to camera vignetting, and changes over time (Faber 1992; 
Baggett & MacKenty 1994). Therefore, galaxies observed 
at different times on different chips and chip positions are 
ideally suited to estimate the uncertainties caused by lack 
of knowledge of the Point Spread Function. Approximate- 
ly 40 % of the total number of objects in the cluster have 
been observed two times, 12 months apart, on different chips 
and / or at different chip positions. Therefore, we are able to 
derive reliable error estimates from a comparison of galaxy 
parameters derived from the two observations. 

A total of 75 objects were observed twice. Fig. 5 shows 
the differences in galaxy parameters derived from the two 
independent sets of observations. Galaxies situated close to 
the edge of a WFC chip are not considered in this compar- 
ison. As expected, the differences in r e and fi e are strongly 
correlated. The slope of the correlation and the spread in r e 
and fi e is somewhat dependent on r e ; the data are grouped 
in three bins of r e to show this behaviour. The slope of the 
correlation is parallel to the line of constant r e I° 8 , which 
enters the Fundamental Plane. There is no systematic off- 
set between the two exposures. With these results, we can 
establish the errors in the determinations of r e and I e . For 
galaxies with r e > 0'.'71 the rms errors in log r e and log I e 
are 0.04 and 0.05, respectively. The rms error for log r e /e ° 8 , 
which enters the Fundamental Plane, is 0.022. 

3.3.5 Comparison with growth curve analysis 

The structural parameters of nearby early-type galaxies are 
usually determined through curve-of-growth fitting. We ver- 
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Figure 5. Comparison of galaxy parameters derived from the 1991 observations (labeled 1) and the 1992 observations (labeled 2). As 
expected, the differences in r e and I e are strongly correlated. The dotted lines are fits to the data. The uncertainty in r e and I e is lower 
for larger (i.e., larger r e ) galaxies. The observed correlations are almost parallel to the lines of constant r e I®' 8 (solid lines), i.e., the 
Fundamental Plane parameter r e I° 8 can be determined to high accuracy. The error is estimated at 0.022 per single observation 



ified that our technique gives similar results. To that end, we 
applied a curve-of-growth analysis to the CLEANed data. Ex- 
cept for the interacting galaxies, the resulting combination 
of r e I° 8 was very similar to that derived from the direct fit. 
The difference was generally below 5 %, except for the SBO 
galaxy in the sample. The median difference in r e I°' & was 2 
%. Hence, the direct fitting method does not introduce large 
systematic effects. 



3.3.6 Recovery of Flux 

We also tested whether the fitting program recovered all 
flux, for stars and galaxies. For stars, we tested this by fitting 
9" x 9" PSFs, generated with Tiny Tim 4.0, to the stars 
used in the calibration of the HST data (cf. Sect. Af). The 
aperture magnitudes were measured through an aperture 
with a diameter of 9" Fig. 6 shows the differences between 
the aperture magnitudes and magnitudes from the fits. The 
mean difference is (_F702PF flt - f 702W 9 ») = 0.02 ± 0.03 
magnitudes. From this, we infer that the results of the fitting 
program and the results from the stellar photometry are 
consistent, and the Tiny Tim construction of the PSF is a 
good approximation of the true PSF. 

Since the wings of the HST PSF extend beyond 3", it 
may be that there are systematic errors in the derived pa- 
rameters I e and/or r e . We tested this by constructing mod- 
el galaxies and fitting the models using PSFs with different 
sizes, 6"x6" and f2"xf2" respectively. The differences in r e 
and I e were smaller than 3 %, whereas the combination of 
r e Ie ' 8 was stable to f %. 

In conclusion, we are confident that the parameters list- 
ed in Table 2 are not systematically different from the true 



galaxy parameters. The typical errors are 0.04 in log r e , 0.05 
in log/e, and 0.022 in logr e / e 8 . 



4 THE FUNDAMENTAL PLANE 

In the late 1980s, it was found that early-type galaxies oc- 
cupy a surface in (r e , I e , <J )-space, now known as the Fun- 
damental Plane (Dressier et al. 1987; Djorgovski & Davis 
1987). The scatter in the edge-on projection of the FP is 
small (< 1 7 % in log r e ), but not entirely explained by mea- 
surement uncertainties (Lucey et al. 1991, JFK95c). The 
source of the intrinsic scatter in the FP is not understood at 
present; the residuals do not correlate with other structural 
parameters like ellipticity or isophotal shape (e.g., JFK95c). 
The form of the Fundamental Plane is given by 

log r e = a log a - /?log I e + 7. (6) 

The precise values of a and fi are not well known; different 
authors find different values, due to differences in sample 
selection and fitting methods. The coefficients depend on 
passband. However, since most data sets are not complete in 
more than one passband, it is often very hard to disentangle 
the effects of differences in passbands and sample selection. 

We therefore use the coefficients measured by JFK95c, 
who studied 225 early-type galaxies in nearby clusters. They 

found «Gunn r = 1-24 ± 0.07 and /fcunn r = 0.82 ± 0.02. 

We use these values throughout the analysis, and we test 
how the conclusions change if we use their coefficients for a 
bluer passband (e.g., «Gunn g = 1.16 ± 0.10 and /?Gunn g = 
0.76 ± 0.04). 

The Fundamental Plane implies that the M/I ratio of 



12 



van Dokkum and Franx 



o 



0.1 



Si 
w 
o 



o - 



-0.1 - 




21 



22 



m 



F702W 



23 

(9") 



24 



Figure 6. Comparison of magnitudes derived with our fitting 
program and magnitudes derived from aperture photometry. The 
fitting program used PSFs with Tiny Tim 4.0 with a size of 9" X 
9" . The aperture magnitudes were derived within 9" . The mean 
difference is 0.02 ± 0.03 magnitudes. We infer that the Tiny Tim 
PSFs are a reasonable approximation of the true PSF. 



galaxies is well-behaved. Under the assumption that galaxies 
are a homologous family, the M/L ratio scales like 



M/L oc (7° 



(7) 



with a scatter of 23 % (Faber et al. 1987, JFK95c). 

The Fundamental Plane is a good tool to study galaxy 
evolution, because the internal scatter is so low. It is there- 
fore sensitive to relatively small changes. In general, evo- 
lution of the Fundamental Plane is expected, because the 
M/L ratios of galaxies are expected to evolve with time. 
This is due to the evolution of the stellar population. Stellar 
populations with a fixed mass are expected to evolve like 



L OC l/(i - iform)" 



(8) 



To first order, k = 1.3 — 0.3x, where x is the slope of the IMF 
(Tinsley 1980). However, k depends also on passband and 
metallicity. As an example, the models of Buzzoni (1989) 
and Worthey (1994) produce values of k between 0.6 and 
0.95 in the V band, for various IMFs and metallicities. 

The evolution of the M/L ratio implies that the zero- 
point of the Fundamental Plane is expected to change. Fur- 
thermore, if the luminosity evolution depends on the mass 
of the galaxy, the coefficients of the Fundamental Plane can 
change (e.g., Renzini & Ciotti 1993). 

The evolution of the M/L ratio is generally measured 
as a function of redshift. Equation (9) implies that the M/L 
evolves in the following way 



where 2f orm is the formation redshift 2f orm of the population 
(Franx 1995, Franx & van Dokkum 1996a, b). This equation 
is valid to first order in z, for go ~ and high values of Z{ ovm . 
The logarithm of the M/L ratio is therefore expected to e- 
volve in a linear way with redshift, and the slope depends 
on a combination of k(IMF), go, and the formation redshift 
2form- In general, the M/L ratio decreases with redshift, be- 
cause the luminosity increases. 



4.1 The Fundamental Plane of CL 0024 

We show the Fundamental Plane for our sample of galaxies 
in CL 0024 in Fig. 7. The galaxies follow a tight relation, sim- 
ilar to nearby cluster galaxies. The scatter is 0.06 in logr e , 
or 15 % in r e . Is is comparable to the observational scatter, 
and similar to the scatter for Coma. Fig. 7a shows the Fun- 
damental Plane along one of its longest axes, whereas Fig. 
7b shows it along a short axis. Similar to nearby clusters, 
it appears tighter in one projection than in the other, but 
the implied scatter is the same for both projections. The 
small scatter in Fundamental Plane implies a scatter in age 
of 20-30 %, depending on the value of k (Renzini & Ciotti, 
1993). 

It is hazardous to determine the coefficients of the Fun- 
damental Plane from the small sample presented here, be- 
cause of systematic errors. Nevertheless, if we use the same 
algorithm as JFK95c, we find coefficients of 



0.97 ±0.09,/? = 0.77 ±0.06 



(10) 



In M/L(z) = In M/L(0) - k (l + go + z' 1 ^) z, 



(9) 



which are not significantly different from the coefficients for 
Gunn g derived by JFK95c. The difference in a is the largest, 
Ad = 0.19 ± 0.14. More, and deeper data are needed to 
measure these coefficients to higher accuracy. If future ob- 
servations confirm that a declines with redshift, then the 
implication is that the M/L ratio of low mass galaxies e- 
volves faster than that of high mass galaxies. This could be 
due to low mass galaxies being younger (as advocated by 
Trager et al. 1993), to differences in the IMF (e.g., Renzini 
and Ciotti 1993), or unknown effects of stellar evolution. 



4.2 The Surface Brightness Test with the 
Fundamental Plane 

We can compare the zeropoints of the Fundamental Plane 
at high redshift and low redshift directly. For the low red- 
shift sample, we use the data of J0rgensen, Franx, & Kjaer- 
gaard (1995a, b) for Coma. Their photometric parameters 
are based on CCD photometry, whereas the velocity disper- 
sions come from a variety of sources. The details can be 
found in the Appendix. The Coma Fundamental Plane has 
a scatter of 16% in the V band. 

Fig. 8 shows the evolution of the zeropoint of the Fun- 
damental Plane expressed in surface brightness. The zero- 
point surface brightness is the surface brightness of a galaxy 
on the Fundamental Plane with an effective radius of 3 
kpc, and a velocity dispersion of 200 km s -1 . The decrease in 
surface brightness is due to the cosmological surface bright- 
ness dimming. Kjaergaard, J0rgensen, & Moles (1993) have 
proposed the Fundamental Plane as a method to verify the 
cosmological surface brightness dimming, but we note that 
the Cosmic Background Radiation provides a much better 
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Figure 7. The Fundamental Plane of CL 0024+ 1654 in the V band, a) An edge-on view of the Fundamental Plane. The two axes are 
two long axes of the Fundamental Plane. The large symbols are the galaxies in CL 0024, the small circles are galaxies in Coma. The open 
square indicates galaxy DGS 162, which is treated in the same way as the other galaxies, but probably is a spiral galaxy. The filled circle 
is the interacting galaxy. The solid line is the fit to the galaxies in Coma by JFK95c. The Fundamental Plane is very similar at z = 0.4. 
b) A similar edge-on Fundamental Plane. The two axes are two of the shortest axes. The Fundamental Plane appears compressed. The 
scatter is more prominent, c) The face-on view of the Fundamental Plane. The galaxies in CL 0024 occupy a smaller region in the Plane, 
which is due to the bias in favor of luminous galaxies. The dashed line is a line of constant luminosity. 
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Figure 8. The evolution of the zeropoint of the Fundamental 
Plane as a function of redshift. The zeropoint is expressed as the 
observed V surface brightness for a galaxy on the Fundamental 
Plane with r e = 3 kpc and <r = 200 km s — 1 , in units of mag 
arcsec - 2 . The surface brightness decreases rapidly with redshift, 
which is due to surface brightness dimming. The drawn line is the 
predicted cosmological effect I e oc (1 + z)~ 4 . The zeropoint for 
CL 0024 lies above this line, which is due to evolutionary bright- 
ening of the stellar population. 

test. If the surface brightness would not evolve like (I + 2) -4 , 
then the CBR would have satisfied a Planck curve for only 
a very short time after the photons were generated. Since it 
satisfies a Planck curve to 0.03 % (Mather et al. 1994), the 
cosmological surface brightness dimming must hold to very 
high accuracy. 

4.3 Evolution of the M/L Ratio Indicated by the 
Fundamental Plane 

Fig. 9 shows the M/L ratios for galaxies in Coma and 
CL0024 as a function of a 0A9 r° e 22 (oc M 02i r-° 02 ). If galax- 
ies evolve in a simple way, without merging or accretion, 
then the parameter a 0A9 r° 22 remains constant, and only 
the M/L ratio will evolve with time. The figure suggests 
that the tilt of the Fundamental Plane may be slightly dif- 
ferent for CL 0024, as found earlier. The reality of the effect 
remains uncertain. 

We can derive the mean evolution in the M/L ratio 
by determining the offsets between both data sets. It turns 
out, however, that a straight fit to the data points produces 
an evolutionary effect which depends strongly on the co- 
efficients of the Fundamental Plane. This is caused by the 
different selection of galaxies. The Coma data are dominated 



by much lower mass galaxies than the CL 0024 data. If the 
coefficients of the Fundamental Plane change, then the Co- 
ma zeropoint will be determined mostly by these low mass 
galaxies, whereas the CL 0024 zeropoint will be determined 
by the high mass galaxies. This difference in the mean mass 
of the galaxies results in a strong correlation between the 
zeropoint difference and the tilt of the lines. The difference 
in (log M) for the two dataset is 0.7, equivalent to a factor 
of 5. The problem can be avoided by giving the low-mass 
galaxies in Coma low weight. All galaxies with masses lower 
than the lowest mass galaxy in CL 0024 were given weight 
0, whereas the other galaxies were given a weight inversely 
proportional to the number of galaxies per mass bin. As a 
result, each mass bin has equal weight, both for the Coma 
sample, and the CL 0024 sample (which was treated in the 
same way). The range in masses is now very similar, and the 
weighted mean masses are the same. 

The fits are shown as lines in Fig. 9a. The mean offset 
between the two datasets is 31% in the M/L ratio. The off- 
set changes very little when different coefficients are used. 
The coefficients for Gunn g give an almost identical offset of 
33%. This can be contrasted to the results when no weight- 
ing is applied: the offsets are 26% and 41% for the two sets of 
coefficients. The uncertainty in the offset due to random er- 
rors is very small, and it is dominated by systematic effects. 
The systematic effects in the structural parameters produce 
an uncertainty of 9 %, whereas the uncertainty due to the 
sample selection is estimated at 7%. The total uncertainty 
is therefore on the order of 12 %. Here we assumed that the 
systematic uncertainties can be added quadratically. 

Figure 9a suggests that the slope for the galaxies in 
CL0024 might be somewhat different from that for Coma. 
This is equivalent to saying that the lower luminosity galax- 
ies have evolved faster in their M/L, which is consistent with 
a change in the coefficient a in the Fundamental Plane. How- 
ever, as we noted earlier, the reality of this effect is not clear, 
and many more data points are needed for verification. 



4.4 Modeling the Evolution of the M/L Ratio 

When we apply the simple formula for a single, co-eval pop- 
ulation to our result, we obtain a constraint of the form 

k(1 + go + l/2 form )0.37 = 0.31 ± 0.12 (11) 

As the models for stellar populations imply 0.6 < k < 0.95, 
the constraint is satisfied for a model with 2f orm = 00 and 
go = 1/2. A model with values of 2f orm = 2 and go = 1/2 
lies 1 a away, for k = 0.6. For high values of k, such models 
predict an evolution of 70%, and that is ruled out firmly. 

These results, however, should be treated with caution. 
As pointed out by Franx (1995), and Franx and van Dokkum 
(1996a, b) the results from the Fundamental Plane may be 
biased if the fraction of early-type galaxies in clusters evolves 
with redshift. In such a case, the set of early-type galaxies at 
low redshift is not the same as the set of early-type galaxies 
at high redshift. Specifically, some of the early-type galax- 
ies at z = may be spirals at z = 0.39, and they should 
be included to derive the proper constraints on the mean 
formation redshift of all early-type galaxies in nearby clus- 
ters. Alternatively, we can try to measure the morphological 
evolution, and correct for it. 
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Figure 9. The evolution of the M/L ratio of early-type galaxies in CL0024. a) The individual measurements of the M/L ratio for 
galaxies in CL 0024 (large symbols), and Coma (small circles). The Fundamental Plane implies that the M/L ratio is proportional to the 
parameter plotted along the horizontal axis. The drawn line is the relation for the CL 0024 sample, the dashed line is the relation for the 
Coma sample. The offset between the Coma sample and the CL 0024 sample is due to the evolution of the M/L ratio. The M/L ratio 
for CL 0024 is 31 % ± 12% lower than that for Coma. The filled circle is the interacting galaxy in CL 0024, the open square is the spiral 
galaxy in CL 0024. There is a suggestion that the M/L ratios of the low luminosity galaxies in C10024 are lower than expected from the 
Fundamental Plane, b) The evolution of the M/L ratio against redshift. The M/L ratio is given for a galaxy on the Fundamental Plane 
with r e = 3 kpc and <r = 200 km s — 1 , in arbitrary units. The hatched area indicates the values allowed by stellar population models 
with £form = 00 and go = 0- The thick lines delineate the area allowed by models with £f rm = 1- The data are marginally consistent 
with such models. 



Another uncertain parameter in the equation is k, the 
luminosity evolution. It depends directly on the IMF, and 
the IMF has not been well determined for early-type galax- 
ies. Hence it could be that the model values of k are incor- 
rect. The measurement of good optical and IR colors may 
help towards a resolution of this issue. The color evolution 
of the galaxies is also linked to the IMF, and the value of k. 
For steep IMF's, k is small, but the color evolution will be 
rapid (e.g., Tinsley 1980). We notice that Stanford, Eisen- 
hard, & Dickinson (1995) found an inconsistency between 
the color evolution of galaxies and the population models. 
The modeling of these observations needs further study. 

Finally, but maybe most fundamentally, we have made 
the assumption that the Fundamental Plane evolves in an 
identical way in Coma and CL 0024. Although there is good 
evidence that the Fundamental Plane is very similar in Co- 
ma, and other rich nearby clusters (e.g., JFK95c), the cluster 
CL 0024 may not be typical of the progenitors of clusters like 
Coma at z = 0.39. Its high richness at z = 0.39 might imply 
that it will evolve into a cluster of much higher richness at 
z = 0. It is therefore valuable to probe the richest clusters 
at z = 0, to test whether the zeropoint of the Fundamental 
Plane is the same. 



5 DISCUSSION AND CONCLUSIONS 

We have measured photometric parameters and velocity dis- 
persions of galaxies in the cluster CL 0024 + 1654 at 2=0.39. 
The Fundamental Plane at this redshift is similar to the 
Fundamental Plane for nearby galaxies. The scatter is smal- 
1, and dominated by measurement errors. Since the scatter 
is expected to increase with redshift, due to the fact that 
fractional age differences between galaxies become larger, 
this is rather surprising. More data are needed to measure 
the shape and scatter of the Fundamental Plane in greater 
detail. The low number of galaxies in our sample makes a 
good determination of these parameters difficult. Further- 
more, the current sample is biased towards the most massive 
galaxies. 

The data demonstrate that massive ellipticals existed 
at 2=0.39 which are very similar to those in nearby clus- 
ters. The formation of these galaxies must have occured at 
significantly higher redshift. This conclusion is in agreement 
with earlier studies on the colors and color magnitude rela- 
tion (e.g., Bower, Lucey, & Ellis 1992). The analysis shows 
a weak evolution of the M/L ratio of the galaxies, on the 
order of 31% ± 12%. Such evolution is expected from the 
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evolution of stellar populations in general. The speed of the 
evolution depends on the IMF, the formation redshift of the 
galaxy, and go- Our low evolution of the M/L ratio is dif- 
ficult to reconcile with a low formation redshift. However, 
large uncertainties in the interpretation remain, relating to 
the IMF, and the progenitor problem: since it is not known 
what the progenitors of present-day ellipticals look like at 
2=0.39, it may be that we are selecting preferentially old 
progenitor galaxies at 2=0.39 and then conclude that they 
are old. Detailed studies of the morphological evolution of 
cluster galaxies is needed to address that issue. 

Our results can be compared to other measurements of 
the luminosity evolution in rich clusters. Aragon-Salamanca 
et al. (1993) derived the evolution of the K band luminosity 
function of rich clusters. They found no significant evolution. 
Their errors are large enough, however, to be consistent with 
the current observations. More such data would be valuable. 
The luminosity function does not necessarily evolve in the 
same fashion as the M/L ratio. If dissipationless merging 
occurs frequently, then its effect on the luminosities is much 
larger than on the M/L ratio. The combined measurement 
of the evolution of the luminosity function, and the M/L 
ratio can constrain such merging. 

We can also compare our results to the evolution of 
the luminosity function of field galaxies. Lilly et al. (1995) 
and Ellis et al. (1995) find little evolution of the luminosity 
function at high luminosities for 2=0.5. This is consisten- 
t with our results. Furthermore, Lilly et al. analyzed their 
red galaxies separately, and found small evolution for the 
red sample as a whole. These results are in qualitative a- 
greement with the current result. A passive evolution of 1 
magnitude, or higher, at 2=0.5 appears to be inconsistent 
with all results. 

It may be more difficult to reconcile the low luminosity 
evolution with the low, or absent color evolution of early- 
type galaxies at 2=0.4 (e.g., Aragon-Salamanca et al. 1993, 
Rakos & Schombert 1995). Stellar population models predict 
rapid color evolution when the luminosity evolution is low, 
and vice versa (Tinsley 1980). 

The datasets therefore agree on a qualitative level, in 
the sense that evolution is small, but more data are needed 
to quantify all evolutionary effects more accurately. Data 
are needed on more galaxies in similar clusters to determine 
the scatter and form of the Fundamental Plane in greater 
detail. Such data can be used to confirm the hint that the low 
mass galaxies have evolved faster than high mass galaxies. 
Furthermore, it will be of interest to study poorer clusters 
at higher redshifts, to test for a dependence of the evolution 
on cluster properties. Finally, the comparison sample at low 
redshift needs to be extended towards the richest clusters. 
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APPENDIX A: TRANSFORMATION FROM 
M FT0 2w MAGNITUDES TO REDSHIFTED V 
MAGNITUDES 

Al Zeropoints 

The HST F702W zeropoint has an uncertainty of order 0.1 
magnitudes (Faber 1992), primarily due to the flatfielding 
problems (cf. Sect. 3.1) and the uncertainty in the contami- 
nation correction. Accurate zeropoints (derived by the Medi- 
um Deep Survey team; Phillips et al. 1994) exist only for the 
F555W and F785LP filters. Hence, we used ground based 
observations of CL 0024 + 1654 to derive the F702W ze- 
ropoint and the transformation to Cousins R and / magni- 
tudes. 

The observations were done on December 12, 1991, on 
the Kitt Peak National Observatory 2.1m telescope, in the 
Cousins R and / bands. The integration times were 1800 s 
in / and 300 s in R. The standard photometric reduction 
was applied, using the IRAF package. The photometric ze- 
ropoints and color terms were calculated with the PHOTCAL 
package. Fluxes of 19 standard stars (Landolt 1992) were 
measured through a 9" diameter aperture. The color terms 
were small, and the published magnitudes could be repro- 
duced to an rms accuracy of 0.02 magnitude. 

The HST data of CL 0024 + 1654 were smoothed to the 



same resolution as the ground based data. The magnitudes 
of the stars and galaxies were measured through a 9" aper- 
ture. We define F702W instrumental magnitudes: 



m,FT02W = 25 — 2.5 log 



ADU 



exptime 
A transformation of the form 

/ = m F 702W - Ol(R — I — (R — /}gal) + /? 



(Al) 



(A2) 



was derived, where (R — /} ga i = 0.82. A least squares fit to 
the data gave the following transformation: 

I = m_F702W-(0.70±0.05)(i?-/-0.82)-(2.92±0.01).(A3) 

The residuals from the fit are shown in Fig. Al. The scatter 
is consistent with the measurement errors. This transfor- 
mation can be compared to the transformations of Harris 
et al. (1991) between the HST filter system and the Kron- 
Cousins filter system, and to the in-flight F702W zeropoint 
derived in Faber (1992) from observations of to Cen. We de- 
rived V — I = 2.13(R — F) from the observed standard stars. 
Hence the equivalent expression to Eq. A3 is given by 



/ = F702W - 0.72(R - I - 0.82) - 2.95. 



(A4) 



The R — I colors of the galaxies fall into a narrow range: 
(R — I) = 0.82 ± 0.06. Therefore, our transformation (E- 
q. A3) is consistent with the transformation from Harris et 
al. (1991) and the zeropoint from Faber (1992) to ~ 0.03 
magnitudes. This difference is well within the uncertainty of 
~ 0.1 magnitudes for the F702W zeropoint (Faber 1992). 

The 9" diameter aperture magnitudes have the advan- 
tage that no aperture correction is needed to compare them 
to the standard star measurements, making them suitable 
for the derivation of Eq. A3. However, they have large er- 
rorbars. 

In the application of eq. A2 to our galaxy photometry, 
we used R — I colors measured through 5" apertures. These 
colors have much smaller errors. The measured values, cor- 
rected for galactic extinction, are listed in Table 2. 



A2 Correction to Redshifted V Magnitude 

The resulting R and / magnitudes can be related to the 
magnitude through a redshifted V filter in the following way. 
Each magnitude is interpreted as a fluxdensity at the central 
frequency v % : m t = c t — 2.5 log F(v t ) The conversion con- 
stants are taken from Frei & Gunn (1994). Their constants 
are based on broad band observations of the spectrophoto- 
metric standards of Oke & Gunn (1983). The uncertainty in 
the constants is determined by the uncertainty in the spec- 
trophotometry. 

Next we assume that the flux density at the redshifted 
V band is related to the flux densities in the R and / band 
by F{v v {z)) = -F(fit) a -F(fj) 1_a . It is straightforward to 
derive the magnitude in the redshifted V band system 



V, 



I + a{R-F) + c I + a{c R - 



-cv + 2.51og(l + z)(A5) 



The last term 2.5 log(l + z) makes the magnitude behave as 
if it is a flux, and not a fluxdensity. This implies that the 
surface brightness decreases as (1 + z)~ l , and not as (1 + 
z)~ 3 as for surface brightness flux densities. The corrected 
magnitude V z behaves therefore as ii'-corrected magnitude. 
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Figure Al. Comparison of ground based / magnitudes with / band magnitudes calculated from HST data using ifi t 



F702W 



0.70(i? — I — 0.82) — 2.92. Filled circles indicate galaxies, open circles indicate stars. In (a), the difference I } 
ground based /. In (b), ig r0 und 
the measurement errors. 



ground — ^fit is plotted against 

Ifil is plotted against R — I color. No systematic trends are apparent. The scatter is consistent with 



We calculated a for galaxy spectral energy distributions 
from Coleman et al. (1980) and Pence (1976), for z = 0.39. 
The value of a is 0.23 ± 0.01 for early-type galaxies and 
early-type spirals. The resulting transformation 



V z = I + 0.23(R - I) + 0.69 



(A6) 



if correct for all galaxy types to 0.01 magnitude. The most 
significant error in the transformation comes from the un- 
certainty in the spectrophotometry. We estimate this uncer- 
tainty at 0.03 from a comparison of synthetic V — I colors of 
stars observed by both Oke (1990) and Hamuy et al. (1994). 



APPENDIX B: THE COMA DATA 

The low redshift comparison sample consists of E and SO 
galaxies in the Coma cluster at z = 0.023. The surface 
photometry is from J0rgensen, Franx & Kjaergaard (1995a) 
[JFK95a]. The data were taken in the Johnson B and the 
Gunn r photometric filters. The effective radii and surface 
brightnesses in V were calculated from the B band photom- 
etry, and B — V colors. The B band is the best band to use 
for this purpose, as the V band is closer to the B band than 
the Gunn r band. Color gradients in B — V are very small, 
and can be ignored. 

The B — V color was calculated from the B — r color. 
The transformation is given by J0rgensen (1994): 



The B — r colors were taken from J0rgensen et al. (1992), 
and seeing corrected as described in JFK95a. 

JFK95a listed the mean surface brightness within an 
effective radius. These are ii'-corrected, and corrected for 
cosmological surface brightness dimming. We removed the 
correction for surface brightness dimming, and converted the 
surface brightnesses to the surface brightnesses at an r e . The 
typical error in the V surface brightness is 0.03 mag. 

The velocity dispersions listed by JFK95b were used. 
These are based on velocity dispersions published by Davies 
et al. (1987), Dressier (1987), Lucey et al. (1991), and 
Guzman et al. (1992). JFK95b corrected these dispersions 
to an aperture size of 3'.'4. 



B -V = 0.673(5 



+ 0.184 



(<r flt = 0.021). 



(A7) 



